Prospects for Detection of Gravitational Waves from Intermediate-Mass-Ratio Inspirals 
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We explore prospects for detecting gravitational waves from stellar-mass compact objects spiraling 
into intermediate mass black holes (BHs) (M ~ 5OM to 35OM ) with ground-based observatories. We 
estimate an event rate for such intermediate-mass-ratio inspirals of < 10-30 yr _1 in Advanced LIGO. We 
show that if the central body is not a BH but its metric is stationary, axisymmetric, reflection symmetric and 
asymptotically flat, then the waves will likely be triperiodic, as for a BH. We suggest that the evolutions 
of the waves' three fundamental frequencies and of the complex amplitudes of their spectral components 
encode (in principle) details of the central body's metric, the energy and angular momentum exchange 
between the central body and the orbit, and the time-evolving orbital elements. We estimate that advanced 
ground-based detectors can constrain central body deviations from a BH with interesting accuracy. 
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First-generation interferometric gravitational-wave 
(GW) detectors, such as such as Laser Interferometer 
Gravitational- Wave Observatory (LIGO) J3] and Virgo J2], 
are searching for GWs at or near their design sensitivities. 
In the next decade, Advanced LIGO (AdvLIGO) H and 
its international partners will increase the volume of the 
Universe searched a thousandfold or more. The most 
promising GW sources for this network are the inspiral 
and coalescence of black hole (BH) and/or neutron star 
(NS) binaries. Current inspiral searches target sources 
with total mass M < 40M Q : NS binaries with masses 
1-3 M Q , BH binaries with masses 3-40 M , and NS-BH 
binaries with components in these mass ranges 

Ultraluminous x-ray observations and simulations of 
globular cluster dynamics suggest the existence of interme- 
diate mass black holes (IMBHs) with masses M ~ 10 2 - 
10 4 M |@]. The GWs from the inspiral of a NS or stellar- 
mass BH into an IMBH with mass M ~ 50-350M Q will 
lie in the frequency band of AdvLIGO. These intermediate- 
mass-ratio inspirals (IMRIs) are analogous to the extreme- 
mass-ratio inspirals (EMRIs) of ~ 1OM objects spiraling 
into ~ 1O 6 M BHs, targeted by the planned LISA obser- 
vatory [7]. We consider NSs and BHs, as less compact ob- 
jects (e.g. white dwarfs) are tidally disrupted at frequencies 
too low to be detectable in AdvLIGO. 

If we consider the possibility that the central body of an 
IMRI (or EMRI) is not a black hole, but some other gen- 
eral relativistic object (e.g. a boson star or a naked singu- 
larity iH), then we can quantify the accuracy with which 
it has the properties predicted for a black hole that: (i) it 
obeys the black-hole no-hair theorem (its spacetime geom- 
etry is the Kerr metric, fully determined by its mass and 
spin), and (ii) its tidal coupling (tide-induced transfer of 
energy and angular momentum between orbit and body) 
agrees with black-hole predictions. Searching for non-BH 
objects may yield an unexpected discovery. 



We report on our initial explorations of the prospects for 
detecting GWs from IMRIs and probing the properties of 
the IMRIs' central bodies. We report on: (i) IMRI event 
rate estimates in AdvLIGO, (ii) estimates of the efficacy of 
GW template families for EVIRI searches, (iii) explorations 
of the character of the IMRI (EMRI) waves if the central 
body is not a black hole, (iv) generalizations of Ryan's the- 
orem concerning the information about the central body 
carried by IMRI and EMRI waves, and (v) estimates of 
the accuracies with which information can be extracted by 
AdvLIGO from IMRIs. 

Event Rates for IMRIs with an IMBH central body. 
We (Mandel et al. UM) estimate that for EVIBH spins 
X = spin angular momentum/Af 2 < 0.3, the distance 
(range) R in Mpc to which a network of three 4 km Ad- 
vLIGO detectors could see IMRIs at a network signal-to- 
noise ratio (SNR) of 8 is 



R 



[1 + (x 2 /2)(M/100M o ) 1 - 5 ] y/m/M Q x 
[800 - 540(M/100M Q ) + 107(M/100M e 



(For IMBHs grown by minor mergers, typical spins will be 
X ~ sJm/M ~ 0.2, with few if any above ~ 0.4.) 

Core-collapsed globular clusters are the most likely lo- 
cations for IMRIs; they may contain an EVIBH and a 
high density of stellar mass BHs and NSs @j. Simula- 
tions show that it is possible to grow IMBHs with masses 
up to M max ~ 350 through a series of mergers 
in the core of a cluster JjJ]. Phinney ifloll estimates an 
upper limit on the EVIRI rate in globular clusters as fol- 
lows: assume each cluster has an IMBH that grows from 
~ 5OM to ~ 35OM by capturing objects of mass m 
in 10 10 years. Core-collapsed clusters have a space den- 
sity of 0.7 Mpc -3 , which gives an estimated IMRI rate of 
~ 0.7 x (300M Q /m) x 10~ 10 Mpc" 3 yr" 1 . This leads to 
a limit of ~ 10 IMRI detections per year in AdvLIGO. 
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A kick velocity Vkick > 50 km/s will eject the merged 
black hole from the cluster, placing an upper limit on m of 
m/M < 0.08 (Vjdck depends on the symmetric mass ratio 
rj = mM/(m + Mf as V kick » I2OOO77V 1 ~ - 
0.937?) km/s 0~D]). Black holes with masses m > 10M Q 
will likely merge with the IMBH or be ejected from the 
core in under 10 10 years. An estimate based on the dy- 
namics of binary hardening via 3-body interactions yields 
a rate of one detection per three years for NS-IMBH in- 
spirals or ten detections per year for BH-IMBH inspi- 
rals lfl2Tl . Optimizing AdvLIGO sensitivity at low frequen- 
cies could improve these rates by a factor of ~ 3. For 
Initial LIGO 0]], rates are much lower due to lower de- 
tector sensitivity and seismic noise below 40 Hz, reducing 
M max to < IOOMq. We estimate an IMRI rate in current 
detectors of < 1/1000 yr^ 1 . 

Search Templates for IMRI Waves with an IMBH 
central body. Matched filter searches require templates of 
sufficient accuracy that the mismatch between template and 
signal does not cause a large loss in event rate. The most 
accurate IMRI templates currently available come from BH 
perturbation theory via numerical solution of the Teukolsky 
equation fH]. Post-Newtonian (PN) template s fill [l5ll and 
PN approximations to Teukolsky waveforms Il6ll are inad- 
equate becuase IMRIs enter the detector frequency band 
when the binary separation is r < 15M and the PN expan- 
sion is poor. 

Inspiral waveforms from black-hole perturbation theory 
are known only to first order in rj plus 0(rj 2 ) in radiation 
reaction. It is important to determine the effect of con- 
servative finite-mass-ratio corrections 0(r? 2 ) but tools to 
study these are not yet in hand. We (Brown fl7ll ) estimate 
these effects by computing the mismatch (for AdvLIGO) 
between restricted PN stationary-phase templates contain- 
ing all known rj terms, and the same templates linearized 
in rj plus 0(r] 2 ) radiation reaction (cut off at the IMRFs 
innermost stable circular orbit); this is the fractional SNR 
loss due to using templates linearized in rj. Mismatches are 
computed at each PN order between 1.0 and 3.5 inclusive. 
For a 1.4 M Q NS-100M o IMBH IMRI, the mismatch is 
< 30% for x < 0.8, and < 15% for X < 0.3. For IMRIs 
with a larger IMBH mass, the mismatch decreases, as ex- 
pected. By allowing the linearized PN waveforms to have 
mass parameters different from those of the nonlinear PN 
waveforms, and minimizing the mismatch over these pa- 
rameters, mismatch falls to less than 10% in all except the 
most rapidly spinning cases I17I1 . Therefore, it is reason- 
able to expect that the Teukolsky waveforms will lose no 
more than 10% of the SNR due to linearization in rj (hence 
no more than a 30% loss of event rate). For detection, it 
will be worthwhile, but not essential, to improve Teukolsky 
waveforms by incorporating nonlinear corrections, but ac- 
curate parameter measurement will require improvements. 

IMRI and EMRI Orbits and Waves; Tri-periodic 
vs. Ergodic. Here we entertain the possibility that the cen- 
tral body is not a black hole. We assume its external space- 



time geometry is stationary, axially and reflection symmet- 
ric and asymptoticaly flat (SARSAF) with metric in the 
formds 2 = — a 2 'dt 2 + w 2 ] {d(j) — ujdt) 2 + gggdO 2 + g rr dr 2 
and all coefficients independent of the Killing time t and 
axial angle <f>. If the spacetime initially is not axisymmet- 
ric, rotation will make it non-stationary; then presumably 
GW emission drives it to stationarity and axisymmetry on 
astrophysically small time-scales. Almost all stationary, 
axially symmetric, self-gravitating objects studied obser- 
vationally or theoretically are reflection symmetric. 

A SARSAF solution to the vacuum Einstein equations 
is determined uniquely by two families of scalar multipole 
moments: mass moments M = M, M 2 (mass quadrupole 
moment), M 4 , . . . ; and current moments Si (spin angu- 
lar momentum), S 3 , S 5 , . . . II 1 811 - For the Kerr metric (de- 
scribing astrophysical black holes), the moments are fully 
determined by the mass M and dimensionless angular mo- 
mentum x = Sx/M 2 via M e + iS e = M l+1 (ixY; this 
is the no-hair theorem. LISA plans to measure as many 
moments as possible, via EMRI waves, and determine the 
accuracy with which each moment satisfies this Kerr for- 
mula; AdvLIGO can do the same for IMRIs. 

For EMRIs and IMRIs, the orbiting object moves along 
an orbit that is nearly a geodesic of the background met- 
ric; gravitational radiation reaction drives it slowly from 
one geodesic to another. If the central body is a Kerr 
BH, then: (i) each geodesic has three isolating inte- 
grals of the motion: energy E, axial angular momen- 
tum L z , and Carter constant Q (and a fourth, "trivial" 
integral, the length of the orbit's tangent vector); (ii) the 
emitted gravitational waves are tri-periodic with W = 

^E« m ^Le' (t!!,+m!! * + " ar) ' ( for inte g er vaIue s of 
k, m, n) il9ll . Here P = +, x is the polarization, and the 

three principal frequencies Qg, fi^, O r , in a precise but 
subtle sense, are associated with the orbital motion in the 
polar (#), azimuthal (eft) and radial (r) directions. The fun- 
damental frequencies and complex amplitudes evolve with 
time as the orbit evolves through a sequence of geodesies. 

If the the Carter constant is lost in SARSAF space- 
times, motion may be ergodic rather than tri-periodic, 
which would make detection of the gravitational waves dif- 
ficult. Gueron and Letelier lEoh have used Poincare maps 
to search for ergodic geodesies in the static (S^ = 0) Erez- 
Rosen metric and we (Gair, et al. 02 ill ) have carried out 
similar studies for a variant of the stationary (Sz ^ 0) 
Manko-Novikov metric l22ll . Both of these metrics have 
arbitrary mass quadrupole moment M 2 , and higher order 
moments fixed by M 2 , Si and M. The Poincare maps 
in these spacetimes reveal that there are geodesies at very 
small radii r ~ few M that appear ergodic, but non at 
large radii. We l2lll found such geodesies only for oblate 
(M 2 < 0) pertubations of spacetimes with sping, but in 
the Erez-Rosen case ll2lll . ergodicity appears only for pro- 
late (M 2 > 0) perturbations. Radiation reaction drives 
the evolution of energy and angular momentum in a way 
that makes it unlikely that the apparently ergodic geodesies 
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could be encountered in the course of an inspiral 112 111 . For 
the apparently non-ergodic (integrable) geodesies, the spa- 
tial coordinates are multi-periodic functions of Killing time 
f to a numerical accuracy of 10~ 7 , and a general argu- 
ment ll23"ll based on the structure of the gravitational prop- 
agator shows that their gravitational waves will have the 
same kind of tri-periodic form as for Kerr BHs. 

There are three possible explanations for the presence of 
large -radius orbits that appear integrable and small-radius 
orbits that appear ergodic in the same spacetime: (i) The 
orbits are actually integrable and actually ergodic, respec- 
tively, (ii) All the orbits are ergodic, but at large radii 
they appear integrable to numerical accuracy because of 
the Kolmogorov-ArnoFd-Moser theorem lE^ll . (iii) All the 
orbits are integrable, but at small radii they are made to 
appear chaotic by some ill-understood numerical instabil- 
ity. It is important to learn which is the case, but for EMRI 
and EVIRI wave observations, apparent integrability (or er- 
godicity) has the same observational implications as actual 
integrability (or ergodicity). 

Information Carried by IMRI and EMRI Waves; 
Generalizing Ryan's Theorem. What information about 
the central body is encoded in the waveforms? We shall as- 
sume the waveforms to be tri-periodic. In principle, a large 
amount of information can be encoded in the time evolu- 
tion of the waves' three fundamental frequencies Clg(t), 
O^t), Q r (t) and their complex amplitudes h Pkrnn (t). It 
has been speculated that these encode, fully and separa- 
bly, the values of all the central body's multipole moments 
{M(,Si} and hence its metric I25I1 . the rates at which 
the orbiting object's tidal pull deposits energy and angular 
momentum into the central body, -Ebody an d Ebody (tidal 
coupling) j26ll , and the orbit's semi-latus rectum pit), ec- 
centricity e(t) and inclination angle i{t) (which carry the 
same information as the isolating integrals) YZ% . That this 



might be so is suggested by a special case that Ryan [25] 
has studied. A trivial extension of Ryan's theorem Il27ll 
leads to the following algorithm for extracting information 
from the waves. Observe the time-evolving modulation 
frequencies as functions of the time-evolving fundamen- 
tal frequency / = f^/V. From these, deduce the func- 
tions Q,A{£l<j,) and thence Q,a{v) for A = 9,r; expand 
in powers of v = (JVffi^) 1 ' 3 ~ (orbital velocity); and 
read out the moments (redundantly) from the two expan- 
sions. Then, knowing the moments and thence the metric, 
use the geodesic equation to deduce p(t) from and 
use wave-generation theory to deduce e(t) and t(t) from 
particular modulation amplitudes, h Pkmn (t). 

We have generalized Ryan's theorem to strongly ellip- 
tical but nearly equatorial orbits (Li 1230), to include tidal 
coupling (Li and Lovelace J27D), and are working on fur- 
ther generalizations. For strongly elliptical but nearly 
equatorial orbits the three fundamental frequencies are in- 
dependent of l at first order. We expand these frequencies 
Qa(M£, Si, e,p) (with A = 9,(j),r) in powers of 1/p, 



with coefficients that depend on e and the moments. Sup- 
pose we observe a series of 2N + 1 values of (£l e , O^, £l r ) 
(for any integer N) during the course of an inspiral. This 
gives us 6N + 3 numbers, from which we can read off 
(via an algorithm based on our expansions of the funda- 
mental frequencies l23ll ): (i) the time evolution of e(t) and 
pit) (2N + 1 values of each), (ii) the lowest + 1 mass 
moments, and (iii) the lowest N current moments. By ob- 
serving the evolving amplitudes of the orbital-precession- 
induced modulations encoded in h Pkrnn , we can recover 
the time evolution of 1. Hence, in principle, we have a full 
description of the spacetime. In practice the methods of 
extracting the information are likely to be quite different 
from these algorithms. 

In the absence of tidal coupling Ryan demonstrated that, 
for a nearly circular, nearly equatorial orbit, the central 
body's moments are encoded not only in the waves' modu- 
lations, but also in the phase evolution of the waves' domi- 
nant harmonic / = f^/V. We have extended this analysis 
to deduce the power being deposited in the central body 
by tidal coupling, -Ebody 1131 ■ We assume the moments 
and metric have been deduced from the precessional mod- 
ulations and then use deviations from the Ryan-theorem 
phase evolution to deduce E hody . Following Ryan, we 

quantify the waves' phase evolution by AN it) = f 2 / f = 
<i(number of wave cycles) /d In /. From this definition 
of AN, we infer the rate of change of orbital energy: 
E orh = idE^/dtt^inl/irAN). All (time-evolving) 
quantities on the right side can be deduced from observa- 
tion plus the geodesic equation (for dE orb / dVt^). From the 
deduced metric and the frequency fit) we can compute the 
power radiated to infinity E^; and thence by energy con- 
servation we can deduce the power being deposited in the 



central body E ho 



dy 



-E, 



orb 



EL, 271. We can also 



infer the angular momentum transferred tidally to the cen- 
tral body, Lbody, via L body = E^y/Q^ (valid for nearly 
circular, nearly equatorial orbits). 

The above argument assumes that we can compute £ m 
without knowing the boundary conditions of the inspiral- 
induced metric perturbation at the central body, since we do 
not know the nature of the central body a priori. For highly 
compact central bodies (those deep inside the perturbing 
field's "effective potential") this is true to high but not com- 
plete accuracy. The effect of boundary conditions at the 
central body on the inspiral phase evolution is communi- 
cated outward to infinity mainly at low frequencies (the 
orbital frequency and its low-order harmonics), and these 
perturbations have great difficulty penetrating through the 
effective potential. If the spacetime metric is Kerr, we have 
shown that the influence of the inner boundary condition 
on the energy radiated to infinity is SE^ ~ v 10 E oo Il27ll — 



five orders smaller in the linear velocity v than the tidal 
coupling -Ebody ~ v 5 E oc ]28|]. Thus, to high accuracy we 



can deduce £ M and thence E bo dy from observations, with 
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out knowing the body's precise nature. 

Measurement Accuracies for AdvLIGO. We have 
estimated how accurately AdvLIGO, via IMRI waves, 
can constrain deviations of the central body's quadrupole 
moment M 2 (Brown [17]) and tidal coupling E'body 
(Fang |29ll ) from those of a Kerr black hole. Absent the 
true waveforms, we used PN waveforms as signals and 
templates. This introduces systematic error, but we be- 
lieve our results are indicative of the accuracies achievable. 
Our source is the circular inspiral of a neutron star into a 
IOOA/q IMBH (under the assumption that radiation reac- 
tion has circularized the orbit (l^A ). The orbit is inclined to 
the hole's equatorial plane, to produce a modulation crucial 
for breaking degeneracy between the IMBH spin \ an d M 2 
and E body . 

To investigate M 2 , we used templates accurate to 3.5PN 
order in phase evolution |fl5ll and 1 .5PN in spin-orbit cou- 
pling II 1411 ad ded the effects of quadrupole-monopole in- 
teraction IBOll to both the phase and the precessional mod- 
ulation and numerically mapped the ambiguity function of 
these signals. For a NS-IMBH IMRI (M 2 = ~X 2 M 3 ) 
with spin \ = 0.8 and SNR ~ 10, we found AdvLIGO 
measurement errors A In M ~ 0.006, A In x ~ 0.02, 
and AlnM 2 ~ 0.6. If the IMBH spin is \ = 0.3, the 
error increases to A In M ~ 0.01, A In x ~ 0.3, and 
A In M 2 ~ 2. The accuracy of measurement depends 
strongly on binary orientation; larger precessional modu- 
lation reduces the errors 03]. 

We model tidal coupling as E hody = eE BH , where 



-E B h is the energy flow into a Kerr black hole 12811 and 
seek to measure deviations parameterized by e. We con- 
structed precessing waveforms l3~ill . with orbital inspiral 
phase given by the 3.5 PN approximation of the Teukolsky 
waveforms iflal . and modulation linearized in inclination 
angle 13211 . We restricted inclination angles to l < 7r/4, 
fixed the direction to the source and the central body's spin 
orientation, and used the Fisher matrix to estimate parame- 
ter measurement accuracies. For a black-hole central body 
with spin x = 0.8 and SNR= 10, we could measure e to 
A In e ~ 1 to 2, increasing to A In e ~ 30 at x = 0.3. 

While these errors are larger than hoped, (i) the 
quadrupole moment \M 2 \ of a boson star with y = 0.3 
is expected to be in the range 15 to 100 % 2 M 3 J8|], so Ad- 
vLIGO could readily identify such a central body, (ii) for 
small spins E Bii — ~\x y5 E oa , and hence for x = 0.3, 
the accuracy of measuring tidal coupling is AE body ~ 
30 x 0.001 Erya, i.e. 3% of the power radiated to infin- 
ity, an interesting accuracy for central bodies with anoma- 
lously large E^ aAy , and (iii) observing an IMRI in each of 
the three AdvLIGO detectors increases the accuracy of pa- 
rameter estimation quoted by a factor of \/3; including ad- 
ditional detectors, e.g. Advanced VIRGO, could improve 
this further. In practice, parameter estimation will be pur- 
sued using Markov Chain Monte Carlo techniques lfl7, 331. 



Our results suggest that AdvLIGO will be able to verify 
with modest but interesting accuracy that an IMRI's central 
body is a black hole, and perform searches for non-Ken - 
central bodies. AdvLIGO's accuracies for probing the cen- 
tral body are far worse than LISA'S (as expected, due to the 
thousand-fold fewer wave cycles), but AdvLIGO is likely 
to be operational some years before LISA. Its studies of 
central bodies will be a valuable precursor to LISA'S EMRI 
science, and might possibly yield a big surprise. 
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